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This work is a follow up of recent investigations, where we study the im- 
plications of a generalized heat kernel expansion, constructed to incorporate 
non-perturbatively the effects of a non- commutative quark mass matrix in a 
fully covariant way at each order of the expansion. As underlying Lagrangian 
we use the Nambu - Jona-Lasinio model of QCD, with SUj(3) and Ua{1) 
breaking, the latter generated by the 't Hooft flavour determinant interaction. 
The associated bosonized Lagrangian is derived in leading stationary phase 
approximation (SPA) and up to second order in the generalized heat kernel 
expansion. Its symmetry breaking pattern is shown to have a complex struc- 
Ph! ture, involving all powers of the mesonic fields allowed by symmetry. The 

considered Lagrangian yields a reliable playground for the study of the im- 
plications of symmetry and vacuum structure on the mesonic spectra, which 
we evaluate for the scalar and pseudoscalar meson nonets and compare with 
^ ■ other approaches and experiment. 
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1. Introduction 

The heat kernel expansion pQ is known as a useful and effective tool to 
study the properties of low-energy QCD 00 El- Depending on the physical 
problem, it can be used either in the form of a derivative expansion [H] , or as 
an inverse mass expansion |Hj. Based on the powerful method of Schwinger - 
DeWitt [7] , it allows for calculations of effective meson Lagrangians directly 
in coordinate space by integrating out the quadratic fluctuations of quark 
fields in presence of a background of classical mesonic fields. The result is cast 
as an asymptotic expansion of the effective action in powers of proper time 
with Seeley - DeWitt coefficients a n , which accumulate the whole dependence 
on the background fields. The remarkable property of the method is that each 
order of the expansion is fully gauge and chiral covariant. 

In the case of massive quantum fields with a degenerate mass matrix M = 
diag(m, m, ...), it is not difficult to derive from the proper time expansion 
an expansion in inverse powers of m 2 , since the mass dependence is easily 
factorized and a subsequent integration over the proper time leads to the 
desired result. The resulting asymptotic coefficients remain unchanged. 

If the mass matrix is however non degenerate M = diag(mi, rri2, ■■■) its 
total factorization is impossible because of the non-commutativity of the ma- 
trix M with the rest of the elliptic operator. It has been shown recently [HI E] 
that masses can be redistributed among the mass-dependent factors by per- 
forming resummations in the series. This leads to new covariant asymptotic 
coefficients. The algorithm for the resummations was derived and the gen- 
eralized heat kernel coefficients b n for the SUf{2) [S] and SUf(3) [5] flavour 
cases were obtained. In ^Uj the relation of the new coefficients with the 
standard ones has been clarified. 

Given the success in the mathematical formulation of the problem, it is 
now a natural step to apply the new asymptotic expansion in the construc- 
tion of effective chiral Lagrangians. This expansion provides a reasonable 
approximation to the physics of massive and heavy quantum fields with a 
non-degenerate mass matrix. This is the case, for instance, of low energy 
QCD. Here a light current quark mass matrix which is non-degenerate is re- 
placed by a non-degenerate mass matrix of heavy constituent quarks through 
the non-perturbative mechanism of spontaneous breakdown of chiral symme- 
try. This area of physics opens a window where our generalized heat kernel 
expansion can be applied. 

Several different approaches based on the standard heat kernel series have 



already been used to study the above mentioned task [21 El El • The main 
difference between them is hidden in the definition of the vacuum state. The 
generalized heat kernel expansion also leads to its own prescription for the 
vacuum. It is clear that the closer one is to the physical vacuum, the more 
realistic the description of the spectrum of the mesonic exitations will be. 
From this point of view we hope that our method is a useful tool for the 
accurate description of the hadronic vacuum state at low energies. 

In the present work we shall choose a well-known quark model |TT] to de- 
scribe the formation of the hadronic vacuum and its mesonic exitations. It is 
an effective microscopic low energy Lagrangian combining the C/z,(3) x Ur(3) 
chiral four-quark Nambu - Jona-Lasinio (NJL) interactions together with the 
't Hooft determinantal six-quark flavour-mixing interaction, responsible for 
£/a(1) breaking [12 a . By including a mass term for the light u, d and strange 
s quarks one can explicitly break the remaining SUl(3) x SUr(3) chiral 
symmetry to the S77/(3) flavour group or its subgroups. This Lagrangian 
has been previously used in [IB*! ITij to calculate the low lying meson mass 
spectrum at leading order 2 . In the recent works ^H] we have analyzed the 
quasi-classical corrections stemming from the 't Hooft interaction, and pre- 
sented a fully analytical solution for the bosonized Lagrangian and effective 
potential. We will use here these results. They represent a necessary step 
for the extension to a larger group of the earlier applications of the method 
in the SU{2) x SU{2) NJL model [HI HE]. 

The paper is organized as follows. In Section 2 we introduce the model 
and present the main results of ^1 needed for the present work. To sum- 
marize, these results are the following. Using path integral methods, the 
bosonization of the fermionic Lagrangian which involves the six quark inter- 
action requires the introduction of two sets of bosonic auxiliary fields, each of 
the scalar and pseudoscalar type, say (s,p) and (c, </>). Then the integration 
over the fermionic fields can be cast in quadratic form, which can be done ex- 
actly. The remaining integrations are over one of the sets of auxiliary bosonic 
variables, (s,p) which are done in the stationary phase approximation. The 
solutions to the stationary path integral equations can be expressed as an 
infinite series in powers of the bosonic scalar and pseudoscalar fields (a, 0), 
with coefficients that are known at all orders. In particular also the symme- 
try breaking piece of the bosonized Lagrangian contains an infinite number 
of terms involving powers of (a, </>) , and are a consequence of the flavour de- 

2 An early approach but without 't Hooft term can be found in ^21 



terminantal interaction. The piece of the bosonized Lagrangian which comes 
from the integration over the fermionic degrees of freedom will be dealt with 
our generalized heat kernel technique in Section 3. In Section 4 we show how 
to deal with the gap equations combined with the requirement of covariance 
of the generalized Seeley - DeWitt coefficients and the symmetry breaking 
pattern of the original Lagrangian which must be not altered. We derive the 
expressions for the masses of the pseudoscalar and scalars in Section 5. In 
Section 6 we present numerical results and conclude with a summary and 
outlook. 

2. The model 

To model low energy QCD, we use the global Ul(3) x Ur(3) chiral symmetric 
four-quark interaction of the NJL-type model 



C 



G 



NJL 



{qKq) + {qiibKq) , (l) 



where A a , a = 0,1,..., 8 are the standard Gell-Mann matrices acting in 
flavour space and normalized by the condition tr(A a Ab) = 25 a b, combined 
with the 't Hooft six-quark flavor determinantal interaction 



£ dct = «(det qP L q + det qP R q) , (2) 

where the matrices Pl,r = (1 -F 7s)/2 are projectors on the left- and right- 
handed quarks. The Lagrangian C& e t lifts the unwanted C/a(1) symmetry 
of £njl for massless quarks, as required by the Ua(X) Adler-Bell-Jackiw 
anomaly of the SUf(3) singlet axial current q^^ri^q in QCD. The total fermionic 
Lagrangian reads 

C = q(iYd„ - rh)q + Ant, (3) 

with the interaction Lagrangian 

Ant = £NJL + £det- (4) 

The quark fields have color (iV c = 3) and flavor (Nf = 3) indices which 
range over the set i — 1,2, 3. The current quark mass, rh, is a diagonal 
matrix with elements diag(m u , rhd,rh s ), which explicitly breaks the global 
chiral SUl(3) x SUr(3) symmetry of the Lagrangian. 

This approach contains several commonly used simplifications which can 
be excluded in a more elaborate consideration. Let us comment first on the 



four-point interaction (JI}. The most general form of this vertex, based on 
phenomenological arguments, needs only to be compatible with the symme- 
try group of low energy QCD and can be chosen to be invariant under the 
SU(3) C x SU L (3) x SU R (3) x U v (l) x U A (1) group. The six-point interaction 
(|2I) corresponds to the N c — > oo limiting case and is modified by the tensor 
term at next to the leading 1/N C order as it follows from the instanton dy- 
namics [TUj . We also assume that all interactions between quarks are taken 
in the long wavelength limit (low momenta) where they are effectively lo- 
cal. The explicit chiral symmetry breaking term, qmq, is standard for QCD. 
There are some doubts in the literature regarding this structure in the con- 
text of the N JL Lagrangian (20] . Further discussion of this point based on 
an instanton approach to the QCD vacuum can be found in pTj . 

In order to access the natural degrees of freedom of low-energy QCD in 
the mesonic sector, we proceed to bosonize the fermionic Lagrangian, by 
introducing in the vacuum persistence amplitude 

Z = fvqVqexp (i fd 4 xC) (5) 

the functional unity J3] 

1 = / Y[^s a Vp a 5{s a - q\ a q)S{p a - qi^Kq) 

HVs a Vp a Va a V<J) a 

^{if***.-™^-^}. (a) 

thus obtaining 

Z = J Yi^o- a V(j) a VqVqexp (i d 4 xC q (q,q,a,(f))) 

x / II ' Ds a'Dpa expii d 4 x£ r (a, (j),s,p)j , (7) 



where 



C q = q{irfdp - a - «750)g, (8) 

G 

£r = — [(s a ) 2 + {p a ) 2 \ + s a (a a - rh a ) + p a <p a 

+ ^A abc s a (s b s c - 3p b p c ) . (9) 



and where the totally symmetric constants A a b c are related to the flavour 
determinant, and equal to 



-™-abc Q. ^ijk^mnl\'^a)im\ / ^b)jn\^c)kl 

= -^d a bc + y — (35ao5bo5 c o — 5 a o5bc — 5bo5 ac — S c oSab 



(10) 



We use the standard definitions for antisymmetric f a b c and symmetric d a bc 
structure constants of U(3) flavour symmetry. One can find, for instance, 
the following useful relations 



feacAbfc + febcAf ac + f e f c A a bc — 0, 

deacAbfc + d e bcAf ac + d e f c A a b c = v6S e0 A ab f. 



;n) 



Here and throughout the paper we use a = cr a X a , and so on for all aux- 
iliary fields, 0, s, p, and use the following representation of the scalar and 
pseudoscalar fields 
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with the following identifications 
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Vn S + vr u 
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K- Rt 



(12) 



V2 J 



77ns - 7T°, <p s = y/2f] s 



On, 



e ns + Oq, &d = e ns — Oq, and a s = ^/2e s for the correctly normalized 
states in the flavour basis (see Eq. (|92)l in Appendix B). Here the subscripts 
ns and s denote non-strange and strange respectively 

For the set of auxiliary mesonic fields s,p the symmetry transformation 
properties are the same as the ones for a, <j) and follow from the chiral trans- 
formations of quark fields 



Sq = i(a + 75/3)?, Sq=-iq(a- f ysP), 



(13) 



where the parameters of the infinitesimal global transformations a and (3 are 
hermitian flavour matrices. One has, for example, 



5s = i[a,s] + {(3,p\, Sp = i[a,p] - {/?, s}. 



(14) 



The symmetry breaking piece of the Lagrangian is contained in C r , since 

6C q = 0, 6C r = 6C SB ^ 0, (15) 

where 

£sb = — -tr(ras) + — -( det(s +ip) + det(s - ip)). (16) 

2 64 v y 

We see that C^b is not invariant under a global chiral transformation due 

to explicit symmetry breaking, governed by the first term, and due to the 't 

Hooft interaction, given by the second term 

5C S B = -tr (ia[m, s] - /3{m,p}) 

+i(3o^(det(s-ip)-det(s + ip)). (17) 

In the following we shall consider the case with diagonal matrix rh where 
rh u = ma 7^ "Ts> i- e - the chiral symmetry is explicitly broken down to 
the vectorial isotopic SUj(2) x U(l)y symmetry. The non-vanishing term 
proportional to k signals Ua(X) breaking leading to the OZI-violating effects 
related to the Adler - Bell - Jackiw anomaly of the SU(3) singlet axial 
current. 

The Fermi fields in Eq.flHJ enter the action bilinearly and the integration 
over them is exact. The result is given in the next section. It is necessary 
to shift the scalar fields in (JJJ), a a {x) — > <J a (x) + m a . It is well known that in 
nature the global chiral symmetry SUl{3) X SUr(3) is spontaneously broken 
down to the Eightfold Way symmetry and the shift takes this into account. 
In the new vacuum state the vacuum expectation values of the shifted fields 
vanish (0|cr tt (x)|0) = 0. The new vacuum is determined by the tadpole mech- 
anism demanding that all tadpole graphs must sum to zero. The constants 
m a denoting the constituent quark masses will be fixed by the gap equations. 

In P^ the lowest order stationary phase approximation (SPA) has been 
used to estimate the leading contribution from the 't Hooft determinant in 
Eq. (J0J) in the functional integrals over s a and p a 

Z[a + m, 4>] = Ml Y\_Vs a Vp a exp (i d A xC r (a + m,4>, s,p)) , (18) 

where M is chosen such that Z[m, 0] = 1. In the SPA the functional integral 
is dominated by the stationary trajectories r" t = (s" t , p" t ), leading to 



Y\ Vs a Vp a exp I i / d 4 x£ r (cr + m, 0, s, p) 



expfi d 4 xC r (r st )j , 



(19) 



where h corrections are neglected. The stationary point, r° t (a, <p\ m), is a 
solution of the equations C r (s, p) = 0: 



Gs a + (a + A) a + — A ahc (s b s c - p b p c ) = 0. 

Gp a + 0a - —KbcSbVc = 0, 

Id 



(20) 



where A a = m a —rh a . This system is well-known from [14. Using expressions 
© and (J2H we obtain 



G r 



6 



v4) 2 + (Pst) 2 ]+^(^ + A )^t + 0aK 



(21) 



One solves Eqs. (|2()jl exactly, looking for solutions s" t and p" t in the form of 
increasing powers in fields a a , 4> a 



h a + h$a b + h { 2 c a h a c + h^.<f> b <f> e + h { ^ cd a b a c a d 



+h {2) 

^"'abcd 



ab u b -r iv abc L 
* + ■■■ 



P 



st 






,(3) 



+ KbaPb<Tc<i>d + k 



abcdMc^d 



+ ... 



(22) 



with coefficients depending on m a and coupling constants. Putting these 
expansions in Eqs. (|2U|1 one obtains a series of selfconsistent equations to 
determine coefficients h a , h ah , /^ 6 and so on. The first three of them are 



3k 

Gh a + A a + — A abc h b h c = 0, 

{ G6ac + T^ Aacbhb ) h ™ = ~ 5a 

G5 ac - —A acb h b j h^J = -S a 



(23) 



All the other equations can be written in terms of the already known coeffi- 
cients, for instance, we have ^EJ 



h 



h 



(i) 

abc 

(3) 
abc 



— h {l) h {1) h {1) A - 

on <*a bb cc abc i 
_^l(2),(2),(1) j, 

. _ i^nn Iwn Ivor -<J-n 



ft 



16 



'aa "'{,6 cc -^abc 1 



(2) _ 
abc 

L abcd 



l^nn K-iI Ai/v^ -*A 



32 

3k 

— / 
16 



aa 6b cc ^aftc ) 

h (l) h (l) h (l)A _ 
"-aa "> bb "'ccd^-abc i 



3 k 



,(2) _ Oft (1) /,(!), (2) _,(2),(3)\ . 
"-afecd _ -i fi ' i aa l/'fto ' i ccd ,l cb n cdb) /i 



abc i 



(24) 



One can see from these equations that the terms quadratic and higher order 
in mesonic fields in Eqs.(|22|) are generated by the 't Hooft interaction and 
will disappear if k — 0. Let us also give the relations following from (|2"3*|) 
which have been used to obtain (J2"3jl 

h b = (Gh a + 2A a )h a 1 b ) = -(3Gh a + 2A a )/$. (25) 

As a result the effective Lagrangian (|21|) can be expanded in powers of 
meson fields. Such an expansion, up to and including the terms which are 
cubic in a a ,(p a , looks like 

1 1 

C r (r st ) = Ka a + -hajja a a b + -h$(j) a <j) b 

+l*a [h^ b a c + (h® + hfl) <p b <p c ] + 0(field 4 ). (26) 

The coefficients h a are determined by couplings G, K and the mean field 
A a . This field has in general only three non-zero components with indices 
a = 0, 3, 8, according to the symmetry breaking pattern. The same is true 
for h a because of the first equation in (|23|) . It means that there is a system 
of only three equations to determine h = h a X a = di&g(h u , h d , h s ), 

Ai + Ght + ^-Y, tijkhjhk = 0. (27) 

Here the totally symmetric coefficients tijk are zero except for the case with 
different values of indices i ^ j ^ k when t uds = 1. The latin indices i,j, k 
mark the flavour states % = u,d,s which are linear combinations of states 
with indices 0, 3 and 8. In Appendix A we collect the matrices which project 
one set to the other and write out exact solutions for Eq.([2HJ). Let us note 
that Eqs.(|27j) must be solved selfconsistently with the gap equations (see 
Eq.([43|) below) to yield the constituent quark masses in leading SPA order. 

3. Heat kernel expansion 

Equation (J26)) contains the piece of the bosonized effective Lagrangian, which 
has no kinetic terms and is obtained in the weak field limit. Now we turn to 
the evaluation of the fermionic functional integral in Eq.(0), which after the 
shift a a (x) — > a a (x) + m a , reads 



Z[Y] = / VqVqexp (i / d 4 x^[i7 M 9 M — (m + a + i^4>)]q\ 



(28) 



where Y collects the background field dependence as indicated below. This 
fermion determinant accounts for the remaining part of the effective La- 
grangian and leads, in general, to nonlocal mesonic vertices with unphysi- 
cal cuts (the quark deconfinement problem). We have resorted here to the 
Schwinger - DeWitt representation for the real part of the corresponding ef- 
fective action, W[y], to obtain in the end the asymptotics for W[Y] in terms 
of local polynomials of background fields and their derivatives given by the 
heat kernel coefficients at coinciding arguments, 

Z[Y]=ex V (W[Y]), 

1 r°° dt 

W[Y] =ln|det£>| = -- / -p(£A 2 )Trexp(-£Z^£ E ), (29) 

2 JO t 

where Tr designates functional trace, the operator D E stands for the euclidean 
Dirac operator in presence of the background fields a, <fi and 

D^D E = m 2 -d 2 + Y, (30) 

with the definition 

Y = il^d^a + 275^0) + a 2 + {m, a} + <p 2 + 175 [a + m, <£}. (31) 

For the regulator p(tA 2 ), needed to keep the integral convergent at t = 0, 
we use two Pauli-Villars subtractions 3 

p(£A 2 ) = 1 - (1 + £A 2 )exp(-£A 2 ), (32) 

where the cut-off A is a free dimensionfull parameter. The regularization 
function p(tA 2 ), being written in terms of a dimensionless variable r = tA 2 , 
fulfills the necessary conditions: p(r) ~ r 2 /2 at r — > and p(r) —* 1 at 
r — > oo. It is important to know to what extent the specific form of this 
function affects our results. It is obvious that the type of used regulator 
does not affect the chiral invariance of the heat kernel expansion, since the 
generalized heat kernel coefficients hi |Hj, which carry the whole symmetry 
properties of the heat kernel expansion, do not depend on it 

my] = -jy^Y, i i-Mk). (33) 

3 A regularization function p must be introduced to define the coincidence limit for the 
Schwinger representation. The regularization of the quark determinant in general should 
be done in accordance with certain requirements (see, for example, the review of R.D. Ball 
in Q])- Some of them are discussed also in |22j . 

10 



Here the expressions for the first four 6j in the case of SU(2)j x U(l)y flavour 
symmetry m u = m^ ^ m s are 

6o = l, 

h = -Y, 



b 2 = — + =^\ 8 Y, 



Y 3 A 2 A 1 

h = ~ + =^X 8 Y - =^X 8 Y 2 - — (dY) 2 , (34) 

3! 6V3 2^3 12 v ; ' v ; 

where we used the definition Ay = m 2 — m 2 . In (|33)l the trace is to be taken 
over colour, flavour and Dirac 4-spinors indices and the regulator-dependent 
integrals Jj are the weighted sums |Sj 

I i =~(2J i (ml) + J i {m*)) (35) 

with 

poo A+ 

J *( m i) = y o ^-P^A 2 )exp(-tm 2 ). (36) 

For the chosen form of the cut-off function we obtain, for instance, 

Urn 2 ) =A 2 -m 2 In (l + ^ J , (37) 

\2 

(38) 

Both of them are divergent in the limiting case A — > 00. 

Thus, the effective Lagrangian depends on the integrals Jj. The more 
terms of the heat kernel series are taken into account, the more the final 
result depends on the form of the cut-off function p{r) and, therefore, the 
more careful one should be choosing a regulator. In the following we restrict 
our study to the two nontrivial terms, 61 and b 2 , in the asymptotic expansion 
of W[Y]. In this case only two integrals, Iq and I\, are involved. If we intro- 
duced in p{r) two independent parameters, instead of one, A, the outcome 
would not depend at all on the form of the regulator, because one can always 
fix these parameters by fixing independently couplings Iq and I± from exper- 
imental data. Actually, we slightly simplified our calculations working with 
only one parameter A, paying for that the price of having some dependence 

11 




on the regularization procedure, which is finally inherited by the constituent 
quark masses. 

The heat kernel series (}3*3*j) defines the asymptotics of the effective action 
for a physical system with the mass matrix m being large compared to the rest 
of the background fields and their derivatives. It corresponds exactly to the 
considered case of low-energy QCD, where the small meson exitations of the 
quark sea take place in the "superconducting" phase with heavy constituent 
quarks. It is interesting to stress that in comparison with the standard Seeley 
- DeWitt coefficients, which transform covariantly with respect to the action 
of the chiral group, our coefficients hi possess more specific transformation 
properties. Indeed, in the broken vacuum state an arbitrary infinitesimal 
variation 5tr(bi), induced by global transformations of the background fields 

5a = i[a, a + m] + {/3, (/>}, 

5(p = i[a,4>] -{f3,a + m}, (39) 

depends on the variation 5Y which is equal to 

5Y=i[a + ^ 5 [3,Y + m 2 }. (40) 

One can see that already the first coefficient b% transforms non-covariantly, 
because m 2 does not commute with a + 7 5 /5 in PH|) . Nevertheless, one can 
prove that £tr(6j) = for all generalized coefficients hi [§]. 

In the present calculations we truncate the heat kernel series at &2- In 
this approximation the effective Lagrangian £ is given by the sum of only 
two local terms £ = £(61) + £(62) + • • • , where 

£(6l) = £ t ad(&l) +£mass(&l), 

£(6 2 ) = £tad(& 2 ) + £kin(& 2 ) + £mass(& 2 ) + £mt(& 2 ). (41) 

Here we distinguish the tadpole terms, £ ta d, from mass terms, £ ma ss, kinetic 
terms, £km, arid interaction terms, £ m t. We have, for instance, 

£tad(&i) = -j-Y[m u {(T u + a d ) + m s a s ], 



12tt 2 



£tad(& 2 ) = - T7 f^A MS [m u (cr„ + a d ) - 2m s a s ]. (42) 



Joined together with the tadpole contribution from Lagrangian (|2T)J) . they 

12 



lead to the gap equations 



Nc 
6n 2 



K + ^m u (3/ - Aush) 



Nc 
6tc 2 



h s + 7r-^m s (3I + 2A US I 1 ] 



0. 



0. 



(43) 



The mass-part of the heat kernel effective Lagrangian contains two con- 
tributions and is given by 



/•(61+62) = NcI ° /V 

mass a 2 \ a 



N c h 



{A us [2V2(3a cr 8 + 



ws ) 



127T 2 
2 1 m 2\J 1 / 2 1 cr_2\_2 1 /-r_2 ™2\„2 



£ + 0?] 



+2(2m 2 + mi )a 2 + (m 2 u + 5m 2 s )a 2 8 + (7m 2 - m 2 s )a 2 

+ {m u + m s ){m u + 2m s )a 2 + (m s - m u ){2m s - m u )(j) 2 f j , (44) 

where we assume that the indices i and / range over the subsets i — 1,2,3 
and / = 4, 5, 6, 7 of the set a = 0, 1, . . . , 8. Thus we have 



0N2 



tf = 27T+7T- + (7T U ) 

a 2 = 2a+a^ + (ag 



0\2 



<p 2 f = 2{K + K- + K°K°), 
a) = 2{K* + K*~ + K*°K*°) 



(45) 



The kinetic term, £^$2), after continuation to Minkowski space, has a 
non-standard factor 



C 



kinl«2j 



N c h 
16tt 2 



tr 



(d,a) 2 + (d, 



It should be rescaled by the redefinition of mesonic fields 

4vr 2 



O n 



9°a, 



$, 9 2 



N c h 



(46) 



(47) 



where the index R stands for the new renormalized fields. 

By virtue of the PCAC hypothesis the coupling g is related to the weak 
decay constants of the pion, /„., or the kaon, fx, 



./; 



m u 



I 



K 



m s + m u 
^9 



(48) 



To see this let us recall Eq.©, where the quarks bilinears q\ a q and g^7sA a g 
have been replaced by the auxiliary fields s a and p a . The SPA aproximation 



13 



used to estimate the path integral over these variables in (J19J1 restricts them 
to the stationary trajectories s® t ,p® t , given by Eq. (j22|) . Thus, we have 

WsKq = Pt, QKq = <* ■ (49) 

The quark operators are finally represented by expansions in increasing pow- 
ers of bosonic fields a a and <fi a . This is a convenient form to establish a 
connection to some current algebra results, such as the PCAC hypothesis or 
the Gell-Mann - Oakes - Renner (GOR) relation J2H|. 
For instance, one easily finds from (|49|) . 

y/2G(l + u a ) V2m u \ g J 

where result ()56|) has been used to obtain the last equality. In exactly the 
same way one derives with the help of Eq. (|57jl 

(A-l-^10) = l f':\f l °l = T^H (==±^) (A'-|^. + |0>. (51) 



Let us assume that (|48|) holds, then these equations coincide with the well 
known PCAC relations. 

One can use the second equation in (J4TJJ) to estimate the quark condensates 
in the vacuum. As far as the isotopic invariance is implemented here we have 

(0|«Ti|0> = (0\dd\0) = ^ , (0|5*|0) = ^ . (52) 

Combining these equations with Eqs.([5l)|). (joTf) and (j4*Sj) one finds the GOR 
relations (up to the last terms in the round brackets, which are proportional 
to the current quark masses and give some model corrections to the leading 
order result) 

1 + ^ J . (53) 

m 2 K f 2 K = -\{m u + m s )(0\uu + ss\0) (l + J^) • (54) 
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4. Mass spectrum 

We proceed now to extract the mass terms for the low-lying pseudoscalar and 
scalar nonets. We discuss first the pseudoscalar spectrum. The quadratic 
terms in the fields from Eq. (|2fij) and Eq. (|ll"|) combine to yield for instance 



•*-"massv^"J 



N 
' (3/o-AJj)- 



12i 21 u Ui 1J 2G(l + u s ) 
rh u <t>1 



2Gm u (l + uo s 



(55) 



To get this result we used the gap equation (|43|) and the stationary phase 
conditions (|27|). Let us also remind that some of our notations and results are 
explained in Appendix A. Finally the pion mass is obtained by introducing 
physical fields (|4*7|) 

2 "• 
2 9 m u fr-r\ 

m - = o^d + u j ' (56) 

In exactly the same way one can obtain the masses of the other members 
of the pseudoscalar nonet 

2 g 2 (m u + m s ) 

mK G(m u + m s )(l + uj u )> {0I> 

1 

m * = ±-\A+B-J(A- Bf + AD 2 ) , (58) 

.2 9 2 



We also have 



m 2 v ,= ?-(A + B + ^{A-B) 2 + 4D 2 . (59) 



fin, fl, 2 UJ a 



A + B= — + — + ri 

m u m s G/i_ 

^ 1 ( hu h s 8u u + u s 
A-B=-\ — - + 



3 \m u m s Gfi. 



Vj2 I hu_ h s 
where [i± — (1 ± u a — 2uj 2 a ). The argument of the square root is 



3 \m u m s G/!_ / 



(A-B) 2 + 4D 2 =(^-^ + -^-) +s(-^-) . (61) 
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It is known that for m u = ma 7^ rn s there is mixing in the 0, 8 channels. 
This part of the Lagrangian has been diagonalized by introducing physical 
fields 1] and 7/ via an orthogonal transformation, as it is discussed in Appendix 
B, with the mixing angle 8 P (in the singlet - octet basis) defined from the 
diagonalization requirement. 

In the limit of vanishing 't Hooft interaction, n = 0, the mixing angle 8 p 
is equal to the ideal one with tan(26 , it j) = 2a/2 and one can conclude that 
V ~ Vns, ff ~ ~f]s- We find in this case 

2 2 9 2 rhu 2 g 2 (m u + m s ) 2 g 2 rh s 

m„ = m n = — , m K = — , m„ = — — . (62) 

Vas Gm u K G(m u + m s ) Vs Gm s 

Using the gap equations one obtains the relations 



m\ - m\ m s m\ - m\ 



2m u (m s - m u ) m u ' 2m u {m s - m % 



1, (63) 



which show the mass splittings within the nonet. 

In the SU(3) limit m u = m^ = m s for non- vanishing k there is no </>o — > 
mixing, since D = 0. One obtains immediately the masses 

,2„ 



™, = m K = m 88 = Gm (64) 



with the singlet-octet mass splitting 



2 2 3g 2uj ff .^ 

m °° - "* = G(l +u)(l- 2u) ' (65) 



where 



nh 1 / / kA u 



is a solution of the stationary phase equation (|2*7|) for the 577(3) case. In the 
chiral limit, rh = 0, the singlet mass moo takes a non-vanishing value. The 
would-be U(l) Goldstone boson receives a mass as a result of the 't Hooft 
interaction. 

We turn now to the scalar sector. The masses of the scalar mesons are 
as follows. For the mesons usually refered to as ao (I G (J PC ) = 1~(0 ++ )) we 
obtain 



hu 1 \ , , 2 2 , a 2 , 2g lo s 
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and for the strange Kq (/(J i 



|(0 + )) we have 



™>h 



G(l 



+ 



U)„ 



m K + Am s m u + 



h u + h s 
m u + m s 



Am s m u 



(68) 



In the 0, 8 channels one must diagonalize the states. Diagonalization pro- 
ceeds as in the pseudoscalar case and the resulting scalar states are denoted 
by e and e' respectively indicating a set of /o (I G (J PC ) = + (0 ++ )) mesons. 
The mixing angle 6 S is defined in the (0, 8) basis. As a result we obtain for 
the corresponding masses 



m c 



m c 



9_ 

2 



A 
A 



B 
B 



(A - B) 2 + AV 2 
{A - B) 2 + AV 2 



(69) 



where 



A 

A 



B=— + — + 



N c h 



m 
h 



m. 



Y2 u $ 



TV" 



3CJ,, 



(m 2 s + ml) + 



u s 



m. 



m, 



V 



\ m„ 771= 



3G[x + 

10s 



u„ 



3Gfi + 



2 + LO s 

G[i + 



(70) 



and 



(A - Bf + AV 2 



(K_ K_ 



w s 



m„ 



m s 



Gfi 



+ 



+ 



uj„ 



Gjj,j 



(71) 



Supposing for a moment that k — 0, we find the mixing angle 6 S to be 
equal ftd, the e-meson is a pure non-strange state, e ns , and the e' is purely 
strange, — e s . The scalar masses become 



m a 
m K, 

m 2 



77* 



77i_ + Am: 



i 



= m K + Am u m s 
m 2 +Am 2 s , 



(72) 
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giving the following mass splittings within the nonet 

m K* ~ m l = 2 ( m s ~ m u)(m s + 2m u ), 

m l B ~ m K* — 2( m s - m u )(2m s + m u ). (73) 

The latter is three times bigger than in the pseudoscalar case 

m l ~ m L = 3 ( m l ~ m U = 6 ( m s - ml). (74) 

Let us consider now the SU(3) limit m u = rrid = m s for k^O. One has 



™lo = ™h = Mi 8 = ml + Ami + „J , 2 x • (75) 



There is no mixing here, since T> — 0, and the singlet state is splitted due to 
the 't Hooft interaction 



Mi - Mi = -— *— . (76) 

00 88 G(l - uo){\ + 2uo) v ; 

Comparing the SU(3) limit of singlet-octet mass splittings in the pseu- 
doscalar, Eq. (jB3j) . and scalar, Eq. ffTo^) . channels, one observes that these ex- 
pressions have opposite signs for the physically reasonable sets of parameters 
(0 < uj < 1/2), where /i_ and /i + are positive. The 't Hooft interaction 
pulls the singlet pseudoscalar state up and the singlet scalar state down with 
respect to the corresponding octet ones. 

To summarize, the pseudoscalar and scalar masses are obtained by means 
of a specific asymptotic expansion 4 of the heat kernel in the framework of 
a simple model for low energy QCD. It can be improved in different ways. 
We have already mentioned some of them in Section 2. Here we also would 
like to point out that in truncating the heat kernel series at second order we 
are neglecting finite size momentum dependent contributions to the one-loop 
fermion determinant that become more important for the heavier particles, 
so that the pole position for extraction of the masses can be modified in a 
sizeable way. However it is well known that the lack of confinement in the 
NJL model introduces serious difficulties with the crossing of non-physical 
thresholds associated with the production of free quark - antiquark pairs, 



4 A summation over all constant meson fields in this series leads to a derivative (long 
wavelength) expansion. 
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which one may encounter by formally continuing the full Euclidean action 
to Minkowski space. These are the main reasons why we decided in this 
simplified version of the model to truncate the series, taking into account 
only the divergent contributions. On one hand, in doing so, we admittedly 
deviate from the original NJL Lagrangian, however in a way which relies 
heavily on its symmetries and asymptotic dynamics, which are fully taken 
into account. On the other hand, this approach gives us, in principle, a 
chance to correct systematically the coefficients Ii of the heat kernel series 
by introducing new parameters in the regularization function p(t, Ai, A2, . . .) 
and fixing them in accordance with phenomenological requirements. This 
procedure, hopefully, can be developed similary to QCD sum rules, like it 
has been done in [23 an d discussed, in particular, in relation with NJL-type 
models in 



5. Numerical results and discussion 

The parameters of the model, m u , rh s , G, k and A are shown in Table 1. 



Table 1 

The main parameters of the model given in the following units: [m] 

[G] = GeV~ 2 , [k] = GeV" 5 , [A] = GeV. 



MeV, 





rh u (m„) 


m s (m a ) 


G 


— K 


A 


a 


4.9 (302) 


167 (519) 


9.3 


0* 


0.95 


b 


2.8 (211) 


85 (356) 


2.8 


157 


1.4 


c 


2.7 (214) 


92 (397) 


3.1 


88 


1.4 


d 


1.2 (171) 


41 (310) 


1.1 


11 


2.3 


e 


0.7 (155) 


24 (296) 


0.6 


1.6 


3.2 


f 


3.2 (227) 


105 (405) 


3.7 


173 


1.3 


g 


4.9 (296) 


161 (493) 


7.6 


664 


0.95 


h 


2.2 (199) 


75 (375) 


2.3 


45 


1.6 


i 


3.6 (242) 


122 (437) 


4.6 


205 


1.2 


J 


3.6 (235) 


109 (382) 


3.7 


422 


1.2 


k 


4.7 (286) 


155 (485) 


7.2 


477 


0.98 


1 


1.5 (179) 


50 (317) 


1.5 


23.4 


2.0 
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Table 2. 

The pseudoscalar nonet parameters in units of MeV (except for the angle 6 p , which 

is given in degrees). 





m n 


m K 


U 


Ik 


m v 


Vflrf 


0v 


a 


138* 


494* 


92* 


125* 


138 


612 


35 


b 


138* 


494* 


92* 


124 


547* 


1504 


2 


c 


138* 


494* 


92* 


131 


526 


958* 


-4 


d 


138* 


494* 


92* 


129 


547* 


1078 


2 


e 


138* 


495* 


92* 


134 


545* 


958* 


2 


f 


137* 


496* 


92* 


128 


532 


1109 


-2 


g 


137* 


496* 


92* 


122* 


507 


1089 


-7 


h 


138* 


495* 


92* 


133 


535 


958* 


-3* 


i 


138* 


495* 


92* 


129 


516 


958* 


-7* 


J 


138* 


494* 


92* 


121* 


547* 


2187 


2 


k 


138* 


494* 


92* 


124* 


497 


958* 


-10 


1 


138* 


494* 


92* 


127* 


547* 


1156 


2 



In Table 2 is the pseudoscalar spectrum, together with the weak decay 
constants /„-, J'k and mixing angle 8 P ; the masses and mixing angle 6 S of the 
scalars are given in Table 3. Inputs are indicated by (*). The latin letter 
labels on the left hand side identify the sets in the tables. 

The following empirical values are taken from [2B|: mj = 139.57018 ± 
0.00035 [MeV], m% = 493.677 ±0.016 [MeV], m v = 547 ±0.12 [MeV], rry = 
957.78 ±0.14 [MeV] for the masses in the low lying pseudoscalar sector. The 
weak decay constants F%* = 130.7 ± 0.1 ± 0.36 [MeV], F« p = 159.8 ± 
1.4 ± 0.44 [MeV] relate to ours through a v2 normalization factor, thus 
fexp _ 9 2.4 MeV and f^ v ^ 113 MeV. 

The scalar masses up to ~ 2 GeV are presently known to be: a (980) = 
984.7±1.2 [MeV], a (1450) = 1474 ±19 [MeV], / (600) = 400-1200 [MeV], 
/o(980) = 980 ±10 [MeV], / (1370) = 1200-1500 [MeV], / (1500) = 1500 ± 
5 [MeV], / (1710) = 1713 ± 6 [MeV], fsT *(1430) = 1412 ± 6 [MeV], where 
the name of the particle is identified with its mass, in order not to clutter 
the notation. In |2Zj there is reported the possibility of existence of a low 
lying strange scalar meson Kq. A broad resonance with mass _K"q(800) = 
797 ± 19 ± 43 [MeV] is observed in J2 
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Table 3. 

The different fits for the masses of the scalar NJL nonet in units of MeV (except of 

the angle 9 S which is given in degrees), as compared with a putative nonet family 

ao(980), -K"o(800), /o(600) and /o(980). The symbols of resonances stand for their 

masses. 





a {a (980)} 


K* {K*(800)} 


e {/o(600)} 


e' {/o(980)} 


S 


a 


620 


933 


620 


1205 


35 


b 


1215 


1164 


346 


1199 


14 


c 


888 


976 


423 


1097 


22 


d 


985* 


968 


249 


1017 


16 


e 


900 


895 


224 


954 


18 


f 


985* 


1050 


441 


1153 


20 


g 


985* 


1150 


601 


1295 


22 


h 


891 


954 


384 


1063 


22 


i 


889 


1021 


489 


1252 


23 


J 


1447 


1346 


399 


1364 


12 


k 


907 


1087 


339 


1248 


24 


1 


1036 


1009 


263 


1053 
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We start the discussion of the scalar and pseudoscalar sectors with the 
following special case shown in set (a). This pattern corresponds to SU(3) 
breaking (rh u ^ m s ) without £/a(1) breaking (k = 0) and has been considered 
in detail in Section 4 (see Eq. fjo^j) for the pseudoscalars and Eq. (|72j) for the 
scalars) . 

The overall description of mass spectra is reasonable, given the simplicity 
of the model. Particular trends are as follows. Fixing m n , ttik, fn, and m v (set 
b) or m v i (set c) to their empirical values, results in reducing the parameter 
K of the 't Hooft interaction by approximately a factor 2 in going from (b) 
to (c) (dropping slightly with increasing value of the cutoff). The masses for 
the scalars and my are highly sensitive to the choice of the r] mass: only 



a 4% reduction of m v value in (b) is needed to get the empirical 



rrir,/ 



(c), 



corresponding however to a 35% drop of the latter with respect to its value 
in (b). Fixing rj to its empirical mass in (b) not only yields a much too heavy 
if, but also too heavy scalars oo, Kq and e' (Table 3). 

Although the order of magnitude for the scalar masses in set (c) is rea- 
sonable, e.g. the mass of ao is obtained within 10% of its experimental value 
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and the Kq mass within 20%, the general trend for a large set of param- 
eters is m ao < tuk* < "V, as opposed to the present empirical evidence 
m K* < fn ao — irif (980)- The latter ordering can be obtained for sufficiently 
low values of k, see set (d), with m ao — m e ' within 2% of the empirical value, 
but at the expense of a very light e and too low values of current and con- 
stituent quark masses. The mass of Kq, being almost degenerate with a , 
remains too large by 20%. 

In set (e) we fix the 5 parameters of the model completely in the pseu- 
doscalar sector, through m^^rriK-, f^^m^^m^. This constrains the n and G 
parameters to comparatively very low values and yields also small quark 
masses; the a and Kq masses are almost degenerate, the Kq mass being 
slightly smaller than the a$ mass. 

In sets (f) and (g) three model parameters are fixed through m^, rax, fn, 
in the pseudoscalar sector and one in the scalar sector m ao , requiring that 
the average value of the 7], rf masses be within 10% of the empirical value. 

In sets (h,i) we fix m^^irtxi fn, m r)' and the mixing angle in the pseu- 
doscalar channels. Results are also quite sensitive to the choice of fx, see for 
instance sets (j) and (k), where the four input values of sets (b) and (c) have 
been kept respectively, fixing the remaining freedom by reducing slightly the 
values of fx- In set (j) a reduction of fx implies an increase in the magni- 
tude of k, increasing the splitting and turning therefore the r( significantly 
heavier [j] remained fixed). The masses of the scalars increase by about 20%, 
as compared to their values in set (b), the lower f a bit less, by 15%. In 
set (k) the reduction of fx implies also an increase in k and therefore in the 
splitting, this time reducing the value of m v (since m v i was kept fixed). The 
splitting in the scalars is also enhanced, the e' is pushed up and e down. The 
masses of ao and Kq increase only slightly. 

In set (1) the input parameters of (b) were kept, but f K chosen larger. 
The parameter k gets reduced and the conclusions are opposite to the ones 
of set (j). 

The values of the mixing angles 6 P and 6 S shown in Tables 2 and 3 are 
consistent with results obtained in |29j in the framework of the linear a model 
with broken U(3) x U(3) symmetry, where 8 p = —5° and 9 S = 21.9°. and with 
the values 8 P ~ 2°, <p s w —14° reported in [30 . The last angle here describes 
the mixing in the flavour basis and corresponds to tp s (see Appendix B) in 
our notations. This agreement is not accidental, since the bosonized NJL 
model is closely related to the linear sigma model [3*H l32j. 
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6. Concluding remarks 

We have analyzed the Nambu - Jona-Lasinio model of QCD in the light of a 
new generalized heat kernel expansion. The result is an effective Lagrangian 
of low energy QCD, incorporating the complete original symmetry pattern, 
but eliminating all non physical thresholds associated with quark- ant iquark 
pair formation due to the lack of confinement of the original Lagrangian. We 
applied the so obtained Lagrangian in the extraction of the low lying spectra 
of pseudoscalars and scalars. The pseudoscalar spectrum turns out to be 
quite satisfactory and we used it partly to fix the main parameters of the 
model. As can be seen from Table 3 the predictions for scalar mesons are 
also not too far from the experimental masses of the lightest known scalars, 
which is remarkable in view of the simplicity of the model. 

There is growing evidence that an isovector ao(980), an isospinor iCj(800), 
as well as two isoscalars / (600) and / (980), are members of the same low- 
lying scalar nonet (321, [3I|, JUS], JS01- There are however different opinions 
about their origin. In our calculation we considered the lightest scalar nonet 
as being qq states. It is in line with ideas presented in |29j . 

The outcome of the model is obtained in the leading order stationary 
phase approximation and can be implemented. There are different sources 
for corrections both at leading order and next to leading order. For instance 
the inclusion of vector and axial-vector mesons can be important for the 
physical picture, because they contribute already at leading order through 
the pseudoscalar - axial-vector and scalar - vector mixings. There are also 
several contributions at next to leading order, e.g. meson loop corrections 
|3TI| and semi-classical corrections to the 't Hooft determinant [16J. As dis- 
cussed in Section 4 of ^0] there are two distinct regimes of chiral symmetry 
breaking, related to small/large six-quark fluctuations. For large fluctua- 
tions the quantum corrections may be numerically relevant. Our aim in the 
present work was to show that the considered new method for the asymptotic 
expansion of the heat kernel, which is in full agreament with all symmetry 
requirements, leads already in its minimal form to realistic results for mass 
spectra. A more detailed description of the scalar nonet in the framework of 
our method, including its decay properties, will be given elsewhere. 

We thank W. Broniowski for discussions and reading our manuscript. 
This work has been supported by grants provided by Fundagao para a Ciencia 
e a Tecnologia, POCTI/35304/FIS/2000 and SFRH/BPD/11579/2002. 
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Appendix A. Consequences of Eq. (|23|) 

The first equation in (|23|) can be written in terms of quark-flavour compo- 
nents hi (see Eq.()27|l). In general the (u, d, s) basis can be transformed to 
the basis (0, 3, 8) by the use of the following matrices u ia and e a i defined as 

m 

y/2 \/2 y/2 \ t / \/2 V3 1 

v^ -V3 , u; ia = -= ^2 -V3 1 I . (77) 



2V3 It t o I V3 i ^ -2 

Here the index a runs a = 0, 3, 8 (for the other values of a the corresponding 
matrix elements are assumed to be zero). We have then for instance h a = 
e a ihi, and hi = uj ia h a . Similar relations can be obtained for A, and A a . 
In accordance with this notation we use, for instance, that h ci = uj ia h^ . 
The following properties of matrices (J77|) are straightforward: u) ia e a j = Sij, 
ZaiUib = S ab , e ai e aj = 5ij/2 and u ia u ie = 25 ae . The coefficients t ijk are 
related to the coefficients A abc by the embedding formula 3uj ia A a j bc e b je C k = 
tijk- The SU(3) matrices A a with index i are defined in a slightly different 
way 2Aj = uji a \ a and A a = 2e a jAj. In this case it follows that, for instance, 
o = a a X a = (Ji\i = diag(a M , a d , a s ), but 2a a A a = o-jAj. 

The solutions of Eq. (|27jl are given in ^B]. One can express all other 
coefficients h a ... in terms of these basic variables. We quote further our result 
for h a b, splitting the range of running indices a, b on three subsets: r,s — 
0, 8, n,m— 1, 2, 3 and /, g — 4, 5, 6, 7, 

, (1,2) _ Sum , (1,2) _ Sfg /yoN 

For the 2x2 matrix with indices 0, 8 we have 

h (1,2) _ - 1 ( 3 =F (4^« - w«) ±>/2(u>u - w 5 ) \ , . 

3G^ ± V ±V2(u u -u.) 3±2(2lu u + lu s )JJ [ > 

with [i± = (1 ± u s — 2uD and 

., : = ^ • (80, 

Quite often the stationary phase equations considered together with the 
gap-equations help us to simplify essentially the results. Here is an useful 
example that shows Eqs. (|27|) and (|43|) at work. 
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Example. Let us consider the expression for the mass of kaons following 
from our mesonic Lagrangian. It is not difficult to obtain that 



2 2 

m K = g 



1 1 fhs_ K_ 

G(l + u} u ) 2 \m s m U/ 



+ {m s -m u f. (81) 



One notices, by using 



^ 2 f--— )=2K-m^), (82) 

which is a direct consequence of the gap equations, that the following relation 
is fulfilled 

-s- — + — + (m, - nitt =g ■ . 83 

2 \m 3 m u ) m u + m s 

Therefore, we obtain 

< = s2 (g(TT^) + ^T^)' (84) 

which can be further reduced to the final result indicated in the Eq. (|57jl . by 
observing that 

This last expression follows immediatly from Eq. (|27|) . 

Appendix B. Diagonalization of the mass matrix and physical states 

To illustrate how the physical fields are chosen in the main part of the text we 
recall here some useful details of the diagonalization procedure and how to 
relate to several different conventions adopted in the literature. Our starting 
point is a quadratic form Q written in the singlet - octet basis (Xq,Xs) 

e = «**>( ££)(£) w 

which can be diagonalized by an orthogonal transformation to the physical 
states (X, X) 

(X\ ( cos 9 sin0 \ (X \ 

\Xj ~~ [ -sin^ costf ) \X 8 )' { ' 
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The angle 6 is extracted from the equation 

2D 

tan2# 



A-B 



(88) 



After some trigonometry the ^-dependence of the diagonalized matrix Q can 
be absorbed in just one term 



Q 



-XX, x) 



( A-B 

1 A+B+ — 

cos 26 









\ 



A + B 



A-B 



'A N 

,a, 



(89) 



cos 



26 I 



It is easy to see that 




A-i 




J -„^t A ~ D \ If a m 2 , ,n 2 


cos 26 "*" \ cos 26 J vv ^ "' 


and therefore 


Q = m 2 x X 2 + m|A 2 , 


2 1 

x 2 

2 1 

m% = — 
x 2 


A + B + sgnl^- ^-)^{A-By + AD* 

\ cos 20 J v 


A + B-sgn(^- j) y/{A-Bf + 4L> 2 

V cos 26 J v 



(90) 



(91) 



Finally, to identify the fields X, X with the physical ones, one should proceed 
as follows. Firstly, find the angle 6 from Eq. (j88|) . choosing, for instance, the 
principal value of arctan26 ) , i.e. — (7r/4) < 6 < (7r/4). Secondly, determine 
the sign of the ratio (A — B)/ cos 2$. Only after having established which 
value of (mx,Trix) is bigger, should one proceed with identification of the 
physical fields, writing down the corresponding rotation ()87j) . 

Alternatively, one can use the nonstrange - strange basis (A ns , X s ), where 



'X w \ J_( V2 1 

.XJ "73 1 -V2 



K X 8 , 



(92) 



Our definition (|87|1 . taken together with Eq. (j92j) . leads to the explicit repre- 
sentation 



X\ I cos ip — sin ip \ /A m 
A"/ I sin'0 cosip ) \X S 



(93) 
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or 

(X, I cosifi sin if) 



(94) 
X) \ — sin-0 costfj J \—X s l 

The angle if) here is equal to if) = 9 + 9a, where 9 ld (9-^ + 9^ = tt/2) is 




determined by the equations sinftd = y2/3, cos^id = l/\/3, and therefore 
ijj = 9 + arctan \[2 ~ 9 + 54.74°. It means that if) is restricted to the range 
9.74° <ip< 99.74°. The angle # = ip - (tt/2) = 9 - 9 id and belongs to the 
interval —80.26° <ip< 9.74°. These two angles correspond to two alternative 
phase conventions for a strange ss-component. 

Here are examples that illustrate the physical interpretation of the given 
formulae, using the results of our calculations obtained in Section 4. 

Example 1. In the case of pseudoscalars with broken SU(3) symmetry but 
without Ua(1) breaking (k = 0) the 0o and 08 components are mixed with 
the angle 9 = 9 ld and (A — B) < 0. Hence, one can conclude from Eq. (|9*T|) 
that the X- state is a heavier one and corresponds to rj'. We have from (|93J) 
r]' = —r] s and rj = ^ ns . However, if Ua(X) symmetry is broken (k ^ 0), one 
has (A — B) > (these are exactly the cases b-1 shown in the Table 2) and 
we must identify the physical fields in opposite order rj' = X, i] = X. 

Example 2. In the case of scalar mesons there is no difference between 
the two patterns k = and k/0. In both cases we have (A — B) < 0, i.e. 
X = e' , X = e. If k = 0, they are pure flavour states X = — e s and X = e ns . 
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